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Abstract 



Discrepancies between Standard Model predictions and experimental mea- 
surements of the fractions R c and Rb of hadronic Z decays to charm and bot- 
tom are investigated. We show that there exists a discrepancy in two com- 
plementary determinations of B{B — > DX). Reducing the branching ratio 
B(D° -» K~tt + ) by ~ 15% from currently accepted values to (3.50 ±0.21)% 
removes the discrepancy. Since B(D° — ► K~it + ) calibrates most charmed 
hadron yields, the reduced value also eliminates the discrepancy between the 
predicted and measured values of R c and mitigates a problem in semileptonic 
B decays. A reduction in B(D° — > K~ir + ) would also mean that roughly 
15% of all D° and D + decays have not been properly taken into account. It is 
shown that if the missing decay modes involve multiple charged particles, they 
would be more likely to pass the requirements for lifetime B tagging at LEP 
and SLC. This would mean that the charm tagging efficiency in Z — > cc has 
been underestimated. As a consequence Rj, would need to be revised down- 
ward, potentially bringing it in line with the Standard Model prediction. 
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I. INTRODUCTION 



Recent discrepancies between theoretical predictions and experimental measurements 
of the fractions R c and Rb of hadronic Z decays to charm and bottom [jlJQ could have 
very serious implications for the Standard Model. It is therefore extremely important to 
determine whether or not there exist possible explanations for these discrepancies which do 
not contradict the current paradigm. To this end we show in Section II that there now exist 
two complementary determinations of B(B — > DX)Q. After making reasonable adjustments 
to charmed baryon yields, we show that the two estimates disagree. We find that a 15% 
reduction in B(D° — > K~tt + ) to (3.50 ± 0.21)% eliminates the problem. We then proceed 
to demonstrate that reducing B(D° — > K~tt + ) also eliminates the discrepancy between the 
predicted and measured values of R c and alleviates a problem in semileptonic B decays. 
One further consequence of the change in B(D° — > K~7T + ) would be that roughly 15% of 
all D° and D + decays have not been properly taken into account. In Section III we show 
that if these missing decay modes involve multiple charged particles, they would be more 
likely to pass the requirements for lifetime B tagging at LEP and SLC In this case, the 
charm tagging efficiency in Z — > cc events would have been underestimated. This would 
necessitate a downward revision in the measured value of Rb which would bring it closer 
to the Standard Model prediction. The effect could in fact be large enough to completely 
eliminate the Rb discrepancy. 



*Throughout this note, CP violation is neglected and for each process its CP-conjugate is implied. 
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II. R c MEASUREMENTS AND B(D° -> K~ir + ) 



The fraction R c of hadronic Z decays to charm, which is predicted to be 0.172 in the 
Standard Model, has recently been measured to be 0.1598 ± 0.0069 |],[|. Similarly, the 
number of charm quarks per B decay (n c ) was historically measured to be smaller than ex- 
pected 0, especially in view of the small measured inclusive semileptonic B decay branching 
ratio. Furthermore, the sum over all branching ratios of exclusive semileptonic B decays 
falls significantly short of the inclusive B(B — > Xiv) measurements |4]]. 

One possible explanation for these discrepancies is that a systematic under-counting of 
charm has taken place. In particular, a common thread in these measurements is a significant 
reliance upon the measured value of B(D° — > K~n + ) to calibrate charm production and 
decay for a wide range of observable decay modes. The CLEO experiment measures || 

B(D° -> K~7r + ) = (3.91 ± 0.19)% , 

and the 1994 Particle Data Group |J cites a world average of 

B(D° -> R-n + ) = (4.01 ± 0.14)% . 

These calibrate not only the D° decay modes, but the D + decay modes as well 0, via the 
ratio 

_ B(D+ -> K-iT + iT + ) 
T+ = B(D° -»■ K-7T+) ' 

The calibration mode for D s , namely B(D S — * 07r), has also recently been tied to B(D° — > 
i^ _ 7r + ) in a model-independent fashion 0. 

As a result of a recent measurement by CLEO of the wrong-charm production in 
flavor-tagged B decays, it is now possible to determine the right-charm branching fraction, 
B(B — > DX), in two complementary ways. As one important consequence, we can treat 
B(D° — > K~n + ) as an unknown which is determined by equating the two results for B(B — ► 
DX). This exercise is carried out in the next section after we address several concerns 
related to charmed baryon yields which result in an overall reduction in the estimate for 
weakly decaying charmed baryon production in B decays. 
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A. Inclusive D Yields in B Decays and B(D° -^K i 



The number of charmed hadrons per B decay is defined as 



n c = Y D + Y Ds + Y baryonc + 2B{B -> (cc)X) 



(2.1) 



where the inclusive production of final states containing an arbitrary charmed hadron T is 
defined by 



Y T = B(B -> TX) + B(B -> TX) 



(2.2) 



The weakly decaying, singly charmed baryon species (A c , 2+'°, fl c ) are collectively de- 
noted by baryon c while (cc) represents charmonia not seen as open charm. Table I sum- 
marizes CLEO measurements with the underlying calibration terms factored out explicitly. 
Note that the inclusive D + yield in B decays involves B{D + — > K~tt + tt + ) which is calibrated 
by D° -> R-7T+ [0 via the ratio, 

_ B(D+ -> K-iT + iT + ) 



B(D° -> iT-vr 

We can thus express Y D + in terms of B(D° — > fC~7r + ) as 

q q% 

y D+ = (0.235 ±0.017) 
= (0.235 ±0.017) 
= (0.238 ±0.029) 
The inclusive D yield in I? decays, 



2.35 ±0.23 . 



(2.3) 



B(D+ -> X-7T + 7T + 

9.3% 
r+ ■ £?(_D° -> K~ix 
3.91% 



5(D C 



7T 



(2.4) 



y D = >>±y D+ 



(2.5) 



can then be expressed in terms of B(D° — > f^~7r + ) as shown in Table I. 

The central values for S c and A c yields which are typically used in the determination of 
n c are both at the 5% level QTOJ . The inclusive A c production in B decays is measured rather 
well, whereas the S c yield has large uncertainty. The CLEO experiment has demonstrated 



that right-sign £ + A c correlations dominate over the wrong-sign i~ A c case ||11|| , (where the 
lepton comes from the semileptonic decay of one B and the A c originates from the other 
B in an T(4S I ) event). As a consequence, inclusive S c production in B decays cannot be 
as large as that of the A c . This is shown in the Appendix where we relate both S c and Q c 
production in B decays to Y\ c and the ratio 

B(B^A C X) , . 

^ S BiB^lx) ' (2 ' 6) 

We neglect b — > u transitions and use the Cabibbo suppression factor, 6 2 = (0.22) 2 , for 

charmed baryon production in b —>■ cus(b —>■ ccd) versus b — > cud'(b — > ccs'). [The prime 

indicates that the corresponding Cabibbo-suppressed mode is included.] The Appendix also 

parametrizes ss fragmentation from the vacuum, and predicts 

-p- = 0.38 ±0.10 , (2.7) 



Yb ar y 0rLc 



B(B — > baryon c X) 



1.41 ±0.12, (2.8) 



1.22 ±0.07, (2.9) 



B(B^baryon c X) = 20±Q1 ^ ^ 

As discussed in the Appendix, S c production in B decay is probably overestimated. Inclusive 
baryon c production thus lies somewhere in the range 

1 < ^ < 1.41 ± 0.12 . (2.11) 

Variation over this range has negligible effect on the value of n c . We therefore use the 
values given in Eqs. (2.7) - (2.10). We also prefer not to use the 1994 PDG value 
of B(A C — ► pK~n + ) = (4.4 ± 0.6)%, because it relies upon a flawed model of baryon 
production in B decays. We instead follow the approach outlined in Ref. and use 
B(A C -> pK~ir + ) = (6.0 ± 1.5)%. Thus n c in Eq. (2.1) can be written : 



n c = (0.883 ±0.038) 
+ (0.042 ±0.008) 



3.91% 



B(D° -> K-71+) 
6% 



+ (0.1211 ±0.0096) 
+ 25(5 -> (cc)X) . 



3.5% 



5(5, 



)7T 



+ 



Inserting the branching fractions in Table II and estimating 13 



B(B -> (cc)X) = 0.026 ± 0.004 



one obtains 



n r = 1.10 ±0.06 



(2.12) 



(2.13) 



(2.14) 



which is below the currently accepted value of 1.18 ± 0.06 [|10| 



Very recently, the CLEO experiment has completed the direct measurement of B(b — ► 
ccs') which allows one to use the following, alternative expression for the number of charm 
quarks per B decay fl~3|j , 



1 — B(b — > no charm) + B(b — > ccs 1 ) . 



(2.15) 



,(-) 



This expression is much less sensitive to either _B(_B— > baryon c X) or B(D° — > K 7r + ). We 
take B(b ^no charm) to be fi~3H , 



B(b -> no charm) = (0.25 ± 0.10) (0.1049 ± 0.0046) 
= 0.026 ±0.010. 



(2.16) 



The inclusive wrong charm B decay branching fraction is expressed as |13|-|l5| 



B(b -> ccs') = B(B -> DX) + B(B -> LTJf) + 



+ B(B ^ baryon c X) + B(B — > (cc)X) 



(2.17) 



From Tables I and III, Eq. (2.17) and the charmed baryon correlations discussed in the 
Appendix, we thus obtain 

3.91% 



B(b^ccs') = (0.085 ±0.025) 

+ (0.100 ±0.012) 
+ B(B^ (cc)X) 



B{D° -> K--K- 
3.5% 

B(D S -> <H 



+ (0.0059 ±0.0031) 



6% 



5(A C -»■ pK-n- 



(2.18) 



Using the absolute charm branching fractions of Table II we obtain 

B(b -> ccs ) = 0.22 ± 0.03 , 



(2.19) 



h r = 1.19 ±0.03 . 



(2.20) 



The quantities n c and n c must be equal. Their difference can be traced to a significant 
discrepancy in two alternative determinations of 5(5 — > DX). On the one hand, one can 
write 

5(5 -> DX) = 1 - 5(5 -> no charm) - 5(5 -> DfX) + 

- 5(5 -> baryon c X) - 5(5 -> (cc)X) . (2.21) 

Inserting the values from Eqs. (2.13) and (2.16) and the current CLEO results, 

3.5% 



B(B -> £>+X) = (0.021 ± 0.010) 

5(5 -> baryon c X) = (0.0365 ± 0.0065) 
into the right hand side of Eq. (2.21) we obtain 



B{D S -> 0vr) 



5(A C ^pK-ir+) 



(2.22) 



(2.23) 



5(5 -> 5X) = (0.89 ±0.02) 



^2.24) 



On the other hand, current CLEO measurements of Y D and r D (see Tables I and III) yield, 

3.91% 



5(5 -> 5X) = (0.798 ± 0.042) 



5(5° -> K-71+) 



(2.25) 



Equating the two determinations of 5(5 — > DX), 



(0.89 ±0.02) = (0.798 ±0.042) 



3.91% 



(2.26) 



5(5° -> X-7T+ 

it follows that either the coefficient (0.798±0.042), or 5(5° -> X-tt+) = (3.91±0.19)%, or 
both, are incorrect. Let us assume for the moment that only 5(5° — > K~n + ) is incorrect. 
We can then solve Eq. (2.26) for 5(5° — > X _ 7r + ) to obtain, 



B(D° -> K-ix + ) = (3.50 ± 0.21)% . (2.27) 

This is considerably smaller than currently accepted values but compatible with the most 
recent measurement from ARGUS ||16|| , 

B(D° -> K~7r + ) = (3.41 ±0.12 ±0.28)%. 

Eq. (2.27), in turn, yields 

B(b -> ccs') = (22.7 ±3.5)% , (2.28) 

n c = fi c = 1.20 ± 0.04 . (2.29) 

Our result must of course be corroborated by additional precision studies. In the mean- 
time we have investigated some consequences of a lower value for B(D° — > K~tt + ). 



B. The Low R c Measurement 

Whereas theory predicts 

T(Z° -> cc) , 

i? c = tt-— ^ — - 1 — r = 0.172 , 2.30 

T(Z° -> hadrons) ' v 7 

experiments yield a combined result which is ~ 2a lower [jl[ 

i? c | e:cp = 0.1598 ±0.0069 . (2.31) 

To analyze this result one must make distinctions among the various contributing mea- 
surements. Those which fully reconstruct a primary D* + are calibrated by B(D° — > K~tt + ). 



These are [17, IS 



R C (DELPHI D*) = 0.148 ± 0.007 ± 0.011 , (2.32) 



R c {OPAL D*) = 0.1555 ± 0.0196 , (2.33) 



with a world-average of 

R C (D*) = 0.150 ± 0.011 . (2.34) 

Unfortunately, the uncertainty in R c measurements due to B(D° — > K~n + ) has not been 
explicitly reported. We therefore conservatively retain the full error on R c to write 

3.84% 



0.172 = (0.150 ±0.011) 



(2.35) 



B(D° -> K-7i+ 

(Note that the calibration factor is different than that used previously because these mea- 
surements have taken the updated PDG value @ B{D° -> K~ir + ) = (3.84 ± 0.13)%.) 
This yields 

B(D° -> R-7i + ) = (3.35 ±0.25)% , (2.36) 

which is near to the value we extracted in Eq. (2.27). 

Note that DELPHI has also measured R c via an inclusive double tag method, where only 
the daughter pion of the D*^ is identified. This method does not involve B(D° — ► K~n + ) 
and the result P1P7|, albeit of limited precision, 

R c (7T + n-) = 0.171+[j-[j^ ± 0.015 (2.37) 

agrees well with the Standard Model. 

Other measurements of R c include a lepton method which has very large systematic 
uncertainties, and measurements from both OPAL and DELPHI that involve direct charm 
counting [p],|19[. The extraction of B(D° — > K~tt + ) from the latter is less straightforward 
since a variety of charmed hadrons are involved. Consider, for instance, the recent OPAL 
result PI, 



i^charm counting) = 0.167 ± OMl(stat) ± 0.011 (sys) ± 0.005 (far) . (2.38) 
OPAL measures 
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R c f(c -> D°) B(D° -> /Tv^) = (0.389 ± 0.037)% , 
i? c /(c -> D + ) B(D + -> ^7r + 7r + ) = (0.358 ± 0.055)% , 
i? c /(c -> L>+) -> 07T+) = (0.056 ± 0.017)% , 

Rc f(c -> A+) 5(A+ -> pK-7T + ) = (0.041 ± 0.020)% . (2.39) 

The fractions are summed by using the updated PDG branching fractions [^] as reference: 





) = (3.84 ±0.13)% 




) = (9.1 ±0.6)% , 




) = (3.5 ± 0.4)% , 


B(A C -»■ JOA' 7T + 


) = (4.4 ± 0.6)% . 



(2.40) 

They assume that the undetected primary S c and Q c production is (15 ±5)% of the primary 
A c production, and thus obtain Eq. (2.38). 

We assume the standard model value R c = 0.172 and again use the more accurate 
estimate for B(A C — > pK~ir + ) of (6.0 ± 1.5)%, rather than (4.4 ± 0.6)%, in order to solve 
for B(D° — > K~n + ). We correlate the inclusive primary production fraction of baryon c to 
that of A c via 

f(c - baryon c ) = f(c - A c ) / (1 - pf , (2.41) 

where p models the production fraction of ss fragmentation relative to // fragmentation 
from the vacuum, (/ = u,d or s) plf . The solution for B(D° — > K~7T + ) is 

(0.00389 ± 0.00037) + a00358±a00055 
B(D° —> K~ir + ) = - - ^± (2 42) 

±jy±^ / p _ (0.00056±0.00017) _ (0.00041±0.00020) ' - 1 
Uc B(D s -K/>n+) (l-p) 2 B(A c ->pK-n+) 

Inserting, 

r+ = 2.35 ±0.23, R c = .172, -> </>tt + ) = (3.5 ± 0.4)% , (2.43) 

and E(A C -> pR-Tx + ) = (6.0 ± 1.5)%, (2.44) 
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we obtain 



B(D° -> K~ir + ) = (3.67 ± 0.36)% . (2.45) 



C. Semileptonic B Decays 

Semileptonic B transitions are among the most intensively studied B decays. They 
consist almost entirely of b — > c£~u transitions, since |V^&/Vd>| ~ 0.1 . Thus a primary lepton 
in B decay is typically accompanied by a charmed hadron. Inclusive semileptonic B decay 
measurements detect the lepton without reconstructing the accompanying charmed hadron. 
As a result, uncertainties from charm are minimal. These measurements also usually involve 
very high statistics and so they are generally very precise 



A variety of semileptonic B decay measurements, where the accompanying charm was 
also seen, have been reported [22]. These include the dominant exclusive B — > D^t~u 



processes, B — > D**{X)t~v transitions, and non-resonant B — > ix X l~ v processes. Com- 
bining all information about semileptonic B decay measurements where the associated charm 
is also seen, one finds a significant shortfall relative to the inclusive measurements 0. De- 
creasing B(D° — > K~7i + ) would alleviate this shortfall, because the semileptonic branching 
fractions with reconstructed charm are inversely proportional to B(D° — > K~n + ) and would 
therefore increase. With some theoretical input we estimate that the value B(D° — > K~tt + ) 
= (2.9 ± 0.4)% eliminates the discrepancy 



D. Summary and Implications 

We have demonstrated that currently accepted values for B(D° —> K~tt + ), namely 

(4.01 ±0.14)% (1994 PDG) 

(3.91 ± 0.19)% {CLEO II) (2.46) 
(3.84 ±0.13)% (1995 PDG update) 
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could be too high. The recent wrong charm B(B — > DX) measurement of CLEO opened up 
a second way to determine right charm B(B — > DX). By equating the two determinations, 
we solved for B(D° — > K~tt + ) to obtain the smaller value, 

B(D° -> K-7V+) = (3.50 ± 0.21)% . (2.47) 

We then demonstrated that reducing the value of B(D° — > K~n + ) enables experimental 
results for R c to agree with theory and diminishes the excess of inclusive semileptonic B 
decays relative to the combined exclusive yields. Table [IV] lists the values of B(D° — > K~tt + ) 
required to eliminate the discrepancy in each of these cases. Combining these values one 
obtains the weighted mean value: 

(B(D° -> K--K+)) = (3.40 ± 0.14)% . 

Additional consequences of a lower value for B(D° — > K~n + ) are discussed in Ref. 
||25|| . We note here that it is possible for a reduction in B(D° — ► K~n + ) to affect the 
discrepancy between theory and measurement for Rt>, the fraction of hadronic Z decays 
to bottom quarks. This connection follows by noting that since B(D° — > K~7T + ) calibrates 
almost all charmed meson branching fractions, a lower value for B(D° — ► fT~7r + ) means that 
a significant fraction of D° and _D + decays have not been observed or properly counted. One 
hypothesis is that these missed decays involve high track multiplicities |26j since such decays 
are more difficult to fully reconstruct due to tracking inefficiencies, particle identification 
errors, combinatoric backgrounds and the presence of undetected neutrals. On the other 
hand, high charged multiplicity decays are more likely to generate a lifetime B tag at LEP 
and SLC since they will more likely yield the high number of significantly displaced tracks 
expected for B decays. We explore this possibility further in the next section. 



III. R b MEASUREMENTS 

Recently the fraction Rb of Z hadronic decays to bb has been measured at LEP [|2"7|-|3~1"|] 
and SLC [|32| using a variety of methods including shape variables, multivariate techniques, 
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high pt leptons, and lifetime tags to distinguish the decays of b quarks from those of lighter 
quarks. While each measurement is consistent, within uncertainties, with the Standard 



Model expectation of R b = 0.2155, they combine to yield R b = 0.2205 ± 0.0016 |33| which 
represents a three standard deviation discrepancy. As seen in Table [V], the highest precision 
contributions to this average are those which use lifetime B tagging. Indeed, the lifetime 
measurements (including the lepton + lifetime result from OPAL) yield a simple weighted 
mean value of R b = 0.2200 ± 0.0017 which dominates the overall result. 

The procedure used for measuring R b is to tag Z — > bb events using any of the above- 
mentioned methods, then subtract backgrounds as estimated from Monte Carlo (MC), and 
estimate the B tag efficiency either by MC or directly from data. Obtaining the B tag 
efficiency from data is more reliable and is possible in all cases where double tagging (tagging 
two B hadron decays in one event) is used. As an illustration of the procedure, if we were to 
ignore backgrounds, the number of tagged hemispheres N t , (where the sphere axis is defined 
by the direction of the highest energy jet) and the number of double tagged events N u would 
be expressed as 

Nt = 2e b R b N z (3.1) 

Nu = C b e 2 b R b N z (3.2) 

For a given B tagging algorithm one counts N t and N tt - Nz is the total number of hadronic 
Z decays. C b is a correlation factor which takes into account the fact that the probability of 
tagging a hemisphere may be correlated with whether or not the other hemisphere is tagged. 
Eqs. (3.1) and (3.2) can be solved for R b and e b : 

2N a 

eb= ON- Z (3 ' 3) 
C b N? . N 

(3 - 4) 

When the B purity of the tagging algorithm is high, so that backgrounds are small, Eqs. 
(3.3) and (3.4) are relatively good approximations. Including backgrounds, one arrives at 
the following generalizations of Eqs. (3.1) and (3.2): 
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N t 

tbRb + e c i? c + e u ds(^ — Rb — R c ) (3-5) 



2N- 



Ntt 2 2 2 

-77- = C b e b R b + C c e c R c + C uds e uds (l — R b — R c ) (3.6) 

As input to these equations one can use either the SM value R c = 0.172 or a value measured 
in a separate analysis, and MC predictions for e c and e u d s ■ A variety of data and MC-based 
studies are performed to estimate the correlation factor C b . C c and C u d s differ negligibly 
from unity. 

We have asked whether there could be an explanation for the experimental result which 
does not contradict the Standard Model. Under the assumption that there could be an 
error in the result obtained for the dominant lifetime measurement technique, it follows 
that this error would have to be common to a number of different measurements and, as 
indicated by equation (3.4), would likely result from an excess of single hemisphere tags 
relative to double-tagged events, or from an overestimate of C b . Due to the many subtle 
differences in the experiments and analyses, it is unlikely that a common error could have 
occurred in the estimates of C b . There could however be a missed or incorrectly estimated 
background that enhances N t and/or suppresses N a . Since the b tagging algorithms used 



by the experiments are of high purity, as shown in Table [VI], we conclude that events with 
two tagged hemispheres should have extremely small non-b backgrounds. By process of 
elimination, one is therefore led to the possibility that some non-b source of single hemisphere 
tags could have been overlooked. From Eqs. (3.4) and (3.5) this would imply that one or 
both of e c and e u d s are wrong. 

The light quark tagging efficiency e u d s can be studied directly with data even for the 
added complication of gluon radiation with splitting to bb or cc |33J. The charmed tagging 
efficiency e c is more difficult to determine since it depends upon the relative populations 
of D°, D + , D s , A c and their decays. We can investigate the sensitivity of R b to e c and 
tuds by plotting the variation in R b as a function of each efficiency with all other quantities 
fixed at their nominal values. For simplicity we set C b = 1 which then results in the following 
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simple expression for R^. 



Rb 



N t 
2N Z 




(3.7) 



Ntt 
N z 



N z 



uds 



For nominal values we choose e u d s = 1.7 x 10 3 and e c = 1.4 x 10 2 which are in the 



R c = 0.172. Finally we use = 0.055 and = 0.012 which are comparable to 

the values in Ref. |HJ] and result in Rb = 0.22 as quoted above for the combined lifetime 
measurements. Figure plots R b function of the change in e c and e u d s • The discrepancy 
between measured and predicted values for Rb is seen to correspond to roughly a 50% increase 
in e u ds from its nominal value or a 20% increase in e c . (The discrepancy is of course also 
removed by smaller but simultaneous upward shifts in both efficiencies.) 

A review of the various measurements of e u d s yields no obvious oversights. Our initial 
concern was that there could be some contribution from gluon radiation followed by splitting 
to bb or cc. However, g — > bb occurs in only about (0.2 - 0.3)% of Z hadronic decays while 
g — > cc occurs in roughly 2.5% of Z hadronic decays P^-|57|. Furthermore, the experiments 
have explicitly studied the effect of g — > bb, cc in Z decays to light quarks and find no 
significant enhancement in e u d s ■ Even under the assumption of a large uncertainty in the 
probability for this phenomenon to occur, it is not possible to achieve anywhere near the 
50% shift required in e uds to explain the Rb discrepancy. 

We next consider e c . We were unable to find explicit consideration of the possible 
effect of gluon radiation and splitting to bb and cc in Z — > cc events. The effect of gluon 
splitting in these events could be more important than in the light quark Z decays since: 
(i) the Rb discrepancy is removed by a smaller change in e c than would be required for 
e u d s and (ii) Z — > ccg with g —>■ cc results in events with typically three heavy flavors in one 
hemisphere. Such hemispheres could be expected to have a higher than average multiplicity 
of significantly displaced charged tracks. The primary charm quark from the Z decay will 
typically remain fairly energetic in spite of having radiated a gluon, and it is then only 
necessary for one of the charm quarks from g — > cc to be energetic in order to have two 



mid-range of the values used by the LEP experiments (see Table |VTp. We assume the value 
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significantly displaced heavy flavor decays similar to that expected for B decay. 

To study this phenomenon we have used the PYTHIA Monte Carlo (version 5.7) with 
JETSET (version 7.4) To cross-check these results we have written a toy MC based 

upon the differential relative rate for the process Z° — ► ccg, g — > cc calculated using an 0(a 2 ) 
cross section. Our calculation is more precise than the leading-log approximation at 0(a 2 ) 
valid only at small angles, because we wanted to also consider large-angle gluon emission. 
We neglected terms proportional to fc 2 /M§, with k 2 denoting the virtuality of the gluon. 
As anticipated, we find that the fraction of those hemispheres with g — > cc that contain a 
large number of displaced charged tracks (e.g. Ndispiaced > 4) increases dramatically. These 
however represent only a small fraction of all hemispheres with g — > cc and upon applying 
the jet probability B tag algorithm |29| we found that, on average, the tagging efficiency 
was not significantly higher than for hemispheres which do not contain g — > cc. 

Finally, as mentioned in the previous section, the accepted value for the branching ratio 
B(D° — ► K~n + ) could be too high. One consequence would be that roughly 15% of D° 
decays would be missing. A similar argument applies to D + decays since the most accurate 
measurement of B(D + — ► K~7t + tt + ) is tied to B(D° — ► K~n + ) |7]]. A number of different 
scenarios could be proposed for the missing decays. We consider here one which could link 
the reduced value for B(D° — > K~tt + ) to the discrepancy between the measured value of 
Rb and the value predicted by the Standard Model. In particular, if the missing decays 
were predominantly multiple charged particle modes, then it can be demonstrated that this 
would increase the value of e c relative to that currently assumed in recent measurements of 
Rb. 

To demonstrate this we have performed a simple Monte Carlo study of the process 
e+e" -> Z ^cc, again using the PYTHIA MC with JETSET. We have chosen the DELPHI 
detector for the purpose of creating a simple model of detector effects such as silicon detector 
acceptance and impact parameter resolution. The results would however be qualitatively 
the same if we were to use OPAL or ALEPH detector acceptances and resolutions. For 
each stable charged particle we have calculated the impact parameter (d) relative to the 
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e + e~ interaction vertex. The calculated impact parameter is smeared according to the 
measured DELPHI resolution function The impact parameter significance is then defined 
as S = d/ad where a<i depends on transverse momentum |39[ . 

For lifetime B tagging we again apply the jet probability algorithm which uses impact 
parameter significance values of charged particle tracks to detect the presence of long-lived 
particles. The resulting hemisphere probability distribution we obtain for e + e~ — > Z — > cc is 



shown in Figure |2| where it is compared to that obtained by ALEPH • The distributions 
are not expected to be in strict agreement since we have not modelled detector efficiencies 
and acceptances in detail. Nevertheless, the agreement is good and should be adequate 
for this discussion. Figure [| plots the hemisphere tagging efficiencies for all D° decays 
containing 4 or more charged particles, and for all D + decays containing 3 or more charged 
particles (normalized to the tagging efficiency for all charm hemispheres), as a function of 
logio(P§ lt ), where Ph is the jet probability obtained for all charged tracks in the hemisphere. 
In measurements of Rf, using this algorithm, a B tag is defined as a hemisphere satisfying 
Ph < Ph 1 * where 10~ 4 < P^* < 10~ 29 for the various experiments. It is clear from Figure 
|3] that the multiple charged track modes have a significantly higher value of e c in this range 
of P H ut values. Despite the presence of fewer charged tracks, the D + efficiency is much higher 
than that of the D° in Figure |||, becaue of its significantly longer lifetime. 

From our previous discussion, we found that the discrepancy is removed by a 20% 
increase in e c . Relative to a nominal value of e c = 0.014, this corresponds to e c = 0.017. Since 
the charm quark hadronizes as a D° or D + roughly 60% and 25% of the time, respectively, 
the ~ 15% of missing decays corresponds to ~ 13% of all charm hemispheres. Let us assume 
that the remaining ~ 87% of charm hemispheres have efficiency e c = 0.014. We will take 
p™ 1 = 10~ 3 - 5 and assume that all missing decays are high multiplicity decays for which 
the efficiencies are given by Figure |3|. This allows us to estimate the maximum impact on 
Rb. From Figure | at Pg 11 = 10~ 3 ' 5 , D + decays to three or more charged particles have 
a tagging efficiency = 5 • e c while D° decays to four or more charged particles have 
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efficiency e° = 1.5 • e c . The adjusted charm efficiency is therefore given by: 

e' c = 0.87 • e c + 0.15 • (0.6 • e° + 0.25 • e+) = 1.20 • e c = 0.017 

Thus, in the extreme where all missing decays are multiple charged particle modes, the Rb 
discrepancy is completely eliminated. We hasten to add that, in reality, charged modes may 
not represent all of the missing decays so that the effect on Rb could be smaller. In fact, 
at the opposite extreme, if the missing decay modes involve few or no charged particles this 
would lead to an increase in the final value of Rb- 

IV. CONCLUSIONS 

We have demonstrated that complementary determinations of B(B — > DX) can be 
reconciled by a downward revision in the branching ratio B(D° — > K~n + ). This revision, 
together with the central role played by B(D° — > K~n + ) in calibrating almost all charmed 
hadron yields, explains the R c discrepancy and diminishes a problem in semileptonic B 
decays. We have speculated that it could also be linked to the R b puzzle. In particular, 
a reduction in B(D° — > K~n + ) would mean that roughly 15% of all D° and D + decays 
have not been properly seen or counted. In the case where all of the missing decays involve 
multiple charged particle final states, we demonstrated that this leads to a higher than 
anticipated lifetime tag efficiency in Z — > cc events which would be adequate to bring the 
measured value of Rb into line with the Standard Model value of 0.2155 . We have also 
explicitly considered g — > cc in Z — > cc events but find that it does not enhance the lifetime 
tagging efficiency. 
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V. APPENDIX: CHARMED BARYON PRODUCTION IN B MESON DECAYS 

Accurate accounting of inclusive charm yields in B decays requires a consistent descrip- 
tion of charmed baryon production, which is lacking in the existing literature. Several 
years ago it was hypothesized that the soft inclusive momentum spectrum of inclusive A c 
production indicates that b — > ccs is the dominant source of A c 's in B decays |40]. The 



hypothesis predicted (i) large wrong-sign £~A C correlations, where the lepton comes from 
the semileptonic decay of one B and the A c from the other B in an T(4S) event; and (ii) 
large S c production in B decays, which at that time had not been observed and was believed 
to be highly suppressed jy]]. Shortly afterwards, CLEO observed the first evidence of H c 
production in B decays (see Table I), but also found that the right-sign £ + A c correlations 



are dominant (see Table III) [O. CLEO measured [p~T 



B(B -> A C X , . 

r A , = — i= ^4 = 0.20 ±0.14 . 5.1 

Because the CLEO measurements of inclusive S c production in B decays involve large 
uncertainties, and their central values appear to us to be too high, this Appendix correlates 
S c and Vl c production in tagged B decays to that of the more accurately measured A c . 
We neglect b — > u transitions and use the Cabibbo suppression factor 6 2 = (0.22) 2 for 
charmed baryon production in b — > cus(b — > ccd) versus b —>■ cud'(b — ► ccs') transitions. 
The parameter p = 0.15 ± 0.05 models the fraction of ss fragmentation relative to // 
fragmentation from the vacuum, where / = u, d or s. (This value for p was chosen to 
demonstrate that even a large value yields a significant reduction in S c production in B 
decays.) 
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We denote by Cud the fraction of 5's which decay to weakly decaying charmed baryons 
which come from b — > cud, and define C 5s , C 5d , Cus analogously. Because our model allows for 
substantial charmless-baryon, charmless-anti-baryon production in 5 decays, C u d is smaller, 
or at most equal to B ad as defined in Ref. |j4"0| . Similar comments can be made for C 5s , C 5d) 
and C us relative to B 5s ,B 5d , and B us - The simplest version of the model predicts 



5(5 - A C X) = (1 - p)(C ud + C- Cd ) (5.2) 

5(5 - A C X) = (1 - p)(C 5s + C- Cd ) (5.3) 

B(B ^ ~ C X) =pC u d + (l- P) Cus + (1 - p) C- cs + p C- Cd (5.4) 

B(B -> E C X) = p (C 5s + C- Cd ) (5.5) 

B(B -> fi c X) = p(C Ss + C 5s ) (5.6) 

5(5 -> H C X) = (5.7) 

The Cabibbo structure 

C S d/(Ccd + Ccs) = Ccd/Ccs' = 6 2 (5.8) 



Cus/ {Cus + C u d) — C us / Cud' — (5.9) 

allows us to express the six observables listed on the left-hand sides of Eqs. (5.2)-(5.7) 
in terms of the two unknowns C u d an d C 5s - The latter are in turn obtained from the two 
measurements involving inclusive A c production in 5 decays, namely Y\ c and r^ cl as follows : 

Cud (1 + A 2 - A 2 r Ac ) 



Y Ac (l- p )(l + A 2 )(l + r Ac ) 

Ccs _ r A c 

C~ d ~ l + A 2 (l-r Ac ) ' 
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(5.10) 



(5.11) 



where 

e 2 



\V CS \ 2 (i-U^ 



(-) (-) _ (-) 

The inclusive H c , f2 c yields in B decays are thus correlated to inclusive A c production, 



B(B -> A C X) 1 



Y Ac l + r Ac 

B(B^A C X) r Ac 



Y Ac l + r A 

B(B — > E C X) C U d f . M . 2 C, 

y =y-\P+0--p)* +T7 

*A C Y Ac < ^ud 

B(B^E C X) Cud Q 



(5.13) 
(5.14) 



( p + (l_ p )Aa + ^L(l-p + pA 2 )} , (5.15) 

^ ^ud > 

p(l + A 2 ), (5.16) 



Y Ac Y Ac c 



ud 



B (B^n c x) _ <W Aa+ cV) j (517) 



5(5 ^ H C X) = . (5.18) 

We have taken p to be a universal quantity and have assumed that the initially pro- 
duced charmed baryon retains its charm [and when applicable, strange] quantum number[s] 
through to its weakly decaying offspring. This is not justified but is conservative in that it 
yields an upper limit for baryon c production in B decays. We typically expect the initially 
produced charmed baryons (via b — > c) to be highly excited, while this is not expected 



of their pair-produced antibaryons (via b — > u or b — > c) [42]. That a sizable fraction of 
these highly excited charmed baryons could break up into a charmed meson, a charmless 
baryon, and additional debris is irrelevant to our discussion which focuses on weakly de- 
caying charmed baryon production in B decays. In contrast, it is important to note that 
H£ — > A C KX could occur significantly [the superscript r denotes excited resonances]. This 
introduces an additional mechanism for A c production in B decays, which may help explain 
the small measured value of r Ac . It also decreases the naive estimate for weakly decaying 
S c production. Because our predictions have not incorporated such effects, they should be 
viewed strictly as upper limits for S c production in B decays. 
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TABLES 



TABLE I. Inclusive Charmed Hadron Production in B Decays as Measured by CLEO 



T 


Y T = B(B -> TX) + B(B -> TX) 


Reference 


D° 


(0.645 ± 0.025) 


3.91% 1 
_B(D'>^K-it+) 




@ 


D+ 


(0.235 ±0.017) [ 


9.3% 

B(D + ->K-TT+TT+ 


)J 


1 


43 


1 


D 


(0.883 ± 0.038) 


3.91% 
_B(D u -^K-ir+) 






D s 


(0.1211 ±0.0096) [s^^j 




1 


44 


1 


A c 


(0.030 ±0.005) 


r 6% 

_B(A c ^ P K-7T+)_ 




1 


45 


1 


'c 


0.020 ± 0.007 


1 


46 


1 




0.028 ±0.012 


1 


46 


1 



TABLE II. Absolute Branching Ratios of Key Charm Decays 



Mode 


BR [in %] 


Reference 


D° -> K--K+ 


3.91 ±0.19 


1 




3.5 ±0.4 


i 


A c — > pK~ir + 


6.0 ± 1.5 


in 



TABLE III. Inclusive Charmed Hadron Production in Tagged B Decays as Measured by CLEO 



Observable 


Value 


Reference 


_ B(B^A C X) 
rAc ~~ B(B^A C X) 


0.20 ±0.14 


[HI 


_ B(B->DX) 
r ° ~ B(B^DX) 


0.107 ±0.034 


1 


_ B(B^D+X) 


0.172 ±0.083 


[47| 



23 



TABLE IV. Extracted Values of B(D° -» K~ir + ). 



Analysis 


B(D° -» ^-vr+) 


B(B DI) 


(3.50 ± 0.21)% 


R C (D ) 


(3.35 ± 0.25)7o 


R c (charm counting) 


(3.67 ± 0.36)% 


Semileptonic BR's 


(2.9 ±0.4)% 


ALL 


(3.40 ±0.14)% 



TABLE V. R b Results as of the Brussels EPS-HEP-95 Conference 



Experiment 


Data Set(s) 


Measurement Type 


R b (R c = 0.172 Fixed) 


ALEPH 


1992 


Lifetime 


0.2192 ± 0.0022 ± 0.0026 


DELPHI 


1992-3 prel. 


Lifetime 


0.2216 ±0.0017 ±0.0027 


DELPHI 


1992-3 prel. 


Mixed 


0.2231 ± 0.0029 ± 0.0035 


DELPHI 


1992-3 prel. 


Multivariate 


0.2186 ± 0.0032 ± 0.0022 


OPAL 


1992-4 prel. 


Lifetime + lepton 


0.2197 ±0.0014 ±0.0022 


ALEPH 


1990-1 


Event Shape 


0.228 ± 0.005 ± 0.005 


SLD 




Lifetime 


0.2171 ± 0.0040 ± 0.0037 


L3 


1991 


Event Shape 


0.222 ± 0.003 ± 0.007 


LEP 




Lepton Fits 


0.2219 ± 0.0039 


LEP+SLD 




Lifetime Fits 


0.2200 ±0.0017 


ALL 






0.2205 ±0.0016 
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TABLE VI. Purity, efficiencies and correlation values. 



Experiment 


B Purity 


e& 




^uds 


c b 


ALEPH 


0.96 


0.26 


1.18x10 


0.88x10 


0.943 


DELPHI 


0.92 


0.21 


1.60xl0~ 2 


2.52xl0~ 3 


0.952 


OPAL 


0.94 


0.23 


1.37xl0~ 2 


l.OlxlO" 3 


1.006 


SLD 


0.94 


0.31 


2.30xl0~ 2 


0.87xl0- 3 


0.995 
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FIGURES 
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FIG. 1. Rf, as a function of the change in the efficiencies e u d s and e c varied separately. 
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FIG. 2. The tag efficiency versus the Log of the probability cut (Logio(Pn)) for Monte Carlo 
e + e~ — > Z — > cc events used in this paper (open circles) as compared with that obtained by ALEPH 
(solid squares) 
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FIG. 3. The ratio of the tagging efficiency for D° hadronic decays containing at least 4 charged 
particles to the decays of all charmed hadrons (black squares) and the similar ratio for D + hadronic 
decay modes containing at least 3 charged particles (open circles) as a function of Logio of the cut 
on hemisphere jet probability Pg 1 *. 
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